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ON STOCHASTIC DIFFERENTIAL EQUATIONS WITH ARBITRARY 
SLOW CONVERGENCE RATES FOR STRONG APPROXIMATION 

ARNULF JENTZEN, THOMAS MULLER-GRONBACH, AND LARISA YAROSLAVTSEVA 


Abstract. In the recent article [Hairer, M., Hutzenthaler, M., & Jentzen, A., Loss of regularity 
for Kolmogorov equations, Ann. Probab. 43 (2015), no. 2, 468-527] it has been shown that there 
exist stochastic differential equations (SDEs) with infinitely often differentiable and globally 
bounded coefficients such that the Euler scheme converges to the solution in the strong sense but 
with no polynomial rate. Hairer et al.’s result naturally leads to the question whether this slow 
convergence phenomenon can be overcome by using a more sophisticated approximation method 
than the simple Euler scheme. In this article we answer this question to the negative. We prove 
that there exist SDEs with infinitely often differentiable and globally bounded coefficients such 
that no approximation method based on finitely many observations of the driving Brownian 
motion converges in absolute mean to the solution with a polynomial rate. Even worse, we 
prove that for every arbitrarily slow convergence speed there exist SDEs with infinitely often 
differentiable and globally bounded coefficients such that no approximation method based on 
finitely many observations of the driving Brownian motion can converge in absolute mean to 
the solution faster than the given speed of convergence. 


1. Introduction 

Recently, it has been shown in Theorem 5.1 in Hairer et al. [D] that there exist stochastic dif¬ 
ferential equations (SDEs) with inhnitely often differentiable and globally bounded coefficients 
such that the Euler scheme converges to the solution but with no polynomial rate, neither in the 
strong sense nor in the numerically weak sense. In particular, Hairer et al.’s work |9] includes 
the following result as a special case. 

Theorem 1 (Slow convergence of the Euler scheme). Let T G (0, oo), d G {4, 5,... }, ^ G M'’*. 
Then there exist infinitely often differentiable and globally bounded functions /r, a: M'’* —)■ M'’* 
such that for every probability space (f2,J^, P), every normal filtration {J^t)t&[Q,T] on (f2,J^, P), 
every standard {iFt)te[o,T]-Brownian motion W: [0,T] x D —)■ M on (D, P), every continuous 

{J^t)t&[o,T]-odapted stochastic process A: [0,T] x D —)■ with Wt G [0,T]: P(A(f) =^-|- 

lo t{X{s)) ds + (t(A(s)) dW{s)) = 1, every sequence of mappings Y"': {0,1,..., n} x D —)■ 

n G N, with Vn G N, A; G {0,1,..., n} : = ^ + [p{Yfi) I + a(Y") (IY((/ + l)T/n) - 

lY(ZT/n))], and every a G (0, oo) we have 

(1) hm {n--E[\\X{T)-Y:\\])=oo. 
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Theorem [^naturally leads to the question whether this slow convergence phenomenon can be 
overcome by using a more sophisticated approximation method than the simple Euler scheme. 
Indeed, the literature on approximation of SDEs contains a number of results on approximation 
schemes that are specihcally designed for non-Lipschitz coefficients and in fact achieve polyno¬ 
mial strong convergence rates for suitable classes of such SDEs (see, e.g., [I^ [TUI EHl UHl El 
EH ElElElEO] for SDEs with monotone coefficients and see, e.g., [HElinilllEllIUlEllllfor 
SDEs with possibly non-monotone coefficients) and one might hope that one of these schemes 
is able to overcome the slow convergence phenomenon stated in Theorem In this article we 
destroy this hope by answering the question posed above to the negative. We prove that there 
exist SDEs with inhnitely often differentiable and globally bounded coefficients such that no 
approximation method based on hnitely many observations of the driving Brownian motion (see 
([^ for details) converges in absolute mean to the solution with a polynomial rate. This fact is 
the subject of the next theorem, which immediately follows from Corollary in Section 

Theorem 2. Let T G (0, cxd), d G {4,5,...}, ^ G Then there exist infinitely often dif¬ 
ferentiable and globally bounded functions p, cr: —)■ such that for every probability space 

(D,J^,P), every normal filtration {iFt)teio,T] on (f2,J^,P), every standard {iFt)te[o,T]-Brownian 
motion W ■. [0, T] x D —)■ M on (D, P), every continuous {iFt)te[o,T]-odapted stochastic process 

X: [0,T] X D ^ wtthVt G [0,T]: F{X{t) = ^ + ia{X{s)) ds + a{X(s)) dW(s)) = 1, 
and every a G (0, cxd) we have 

(2) lim (^n" ■ inf inf E 

Tl —^CXD \ Sp . ,5^1 G [OjX'] U'. 

measurable 

Even worse, our next result states that for every arbitrarily slow convergence speed there 
exist SDEs with inhnitely often differentiable and globally bounded coefficients such that no 
approximation method that uses hnitely many observations and, additionally, starting from 
some positive time, the whole path of the driving Brownian motion, can converge in absolute 
mean to the solution faster than the given speed of convergence. 


\\X{T)-u{W{si),...,W{sn))\\ ) =oo. 


Theorem 3. Let T G (0, cxd), d G {4, 5,... G and let {an)n&n C (0, cxd) and (dn)neN C 
(0, oo) be sequences of strictly positive reals such that lim„_,.oo Un = hm„_,.oo dn = 0. Then 
there exist infinitely often differentiable and globally bounded functions /i, a: —)■ such 

that for every probability space (D,J^, P), every normal filtration (J^i)tg[o,T] on (D,J^, P), ev¬ 
ery standard {iFt)t£io,T]-Brownian motion IE: [0,T] x D —)■ M on (f2,J^,P), every continu¬ 
ous {iFt)t£[o,T]-odapted stochastic process X : [0,T] x D —)■ with 'it G [0,T]: P(X(f) = 


^ + /o ds + /g cr(X(s)) dW{s)) = 1, and every n gN we have 


(3) inf inf E 

si,.--)-5n£[0,T] u: R^xC([(5^,T])— 


|X(T) - n(lE(si),..., lE(s„), (lE(s)),g[,„,T: 


> 


measurable 


Theorem is an immediate consequence of Corollary in Section together with an appro¬ 
priate scaling argument. Roughly speaking, such SDEs can not be solved approximately in the 
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strong sense in a reasonable compntational time as long as approximation methods based on 
finitely many evalnations of the driving Brownian motion are nsed. In Section we illustrate 
Theorems and 1^ by a numerical example. 

Next we point out that our results do neither cover the class of strong approximation algo¬ 
rithms that may use finitely many arbitrary linear functionals of the driving Brownian motion 
nor cover strong approximation algorithms that may choose the number as well as the location 
of the evaluation nodes for the driving Brownian motion in a path dependent way. Both issues 
will be the subject of future research. 

We add that for strong approximation of SDEs with globally Lipschitz coefficients there is a 
multitude of results on lower error bounds already available in the literature; see, e.g., mini 
I2ni ED ED ED E5]) the references therein. We also add that Theorem 2.4 in Gyongy [7] 
establishes, as a special case, the almost sure convergence rate 1 / 2 — for the Euler scheme and 
SDEs with globally bounded and infinitely often differentiable coefficients. In particular, we 
note that there exist SDEs with globally bounded and infinitely often differentiable coefficients 
which, roughly speaking, can not be solved approximatively in the strong sense in a reaonsable 
computational time (according to Theorem above) but might be solveable, approximatively, 
in the almost sure sense in a reasonable computational time (according to Theorem 2.4 in 
Gyongy [7]). 


2. Notation 

Throughout this article the following notation is used. For a set A, a vector space V, a 
set B ^ V, and a function f: A ^ B we put supp(/) = {x & A\ f[x) 7 ^ 0}. Moreover, for 
a natural number d G N and a vector n G we denote by ||n||]Rd the Euclidean norm of 
V G Furthermore, for a real number x G M we put [xj = max(Z fl (—oo,x]) and [x] = 
min(Z n [x, C)o)). 

3. A FAMILY OF STOCHASTIC DIFFERENTIAL EQUATIONS WITH SMOOTH AND GLOBALLY 

BOUNDED COEFFICIENTS 

Throughout this article we study SDEs provided by the following setting. 

Let T G (0, 00 ), let (D, P) be a probability space with a normal filtration (J^t)tg[o,'r], and 
let W ■. [0, T] X D — )■ M be a standard (J^i)ig[o,T]-Brownian motion on (D, P). Let Ti, r 2 , Ta G M 

satisfy 0 < Ti < r 2 < ra < T and let f,g,h G ^“(M, M) be globally bounded and satisfy 
supp(/) C (- 00 ,Ti], inG 6 [o,n/ 2 ] |/'(s)| > 0, supp(^) C [r 2 ,ra], 4|^(s)|^ds > 0, supp(h) C 
[ra, cx)), and h{s) ds 7 ^ 0 . 

For every G C'°°(]R, (0, 00 )) let /i^: —)■ and a: —?• be the functions such that 

for all X = (xi,..., X4) G we have 

(4) /i'^(x) = ( 1 , 0, 0, h(xi) • cos(x 2 'ipix^))) and a{x) = (O, f{xi),g{xi), O) 
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and let = (Xf ,..., Xf): [0,T] x n ^ be an (j;)te [o,r]-adapted continuous stochastic 
processes with the property that for all t G [0, T] it holds P-a.s. that X^{t) = p'^(X’^(s)) ds + 

f‘^(x%,))dw(sy 

Remark 1. Note that for all V’ G (0, oo)) we have that and a are inhnitely often 

differentiable and globally bounded. 


Remark 2. Note that for all ip G ^“(M, (0,oo)), t G [0,T] it holds P-a.s. that 


A-*(() = (, X^(t) = 


h/ 


f{s) dW{s), 


'0 


(5) 


xtw = lh.T|(«) 


r*min{t,r3} 


g{s) dW{s), 


' min{t,r2} 

Xt{t) = l[r 3 ,T](t) ■COs{Xt{Ti)lp{X^{T^))) ■ I h{s) ds. 


f 7-3 


Example 1. Let Ci, C 2 , C 3 G M and let f,g,h: M —)• M be the functions such that for all x G 
we have 


( 6 ) 


f(x) = l(_,j„ 7 ,)(a:) ■ exppi H-^—), 

V X — Tt / 

g{x) = l(r 2 ,T 3 ){x) ■exp(^C 2 + 
h{x) = 1 ( 7 - 3 , 00 ) (t) ■ exp(^C 3 


X — Ti- 

1 1 
- + - 


T 2 - X X -T 3 
1 


T3 -X. 

Then /, g, h satisfy the conditions stated above, that is, /, g, h are inhnitely often differentiable 
and globally bounded and f^g^h satisfy supp(/) C (—cx),ri], inf^gjorpa] |/'(s)| > 0 , supp(^) C 


[t 2 , r 3 ], /k |^(s)|^ ds > 0, supp(h) C [rs, cx)), and C Ks) ds ^ 0. 


4. Lower error bounds for general strong approximations 


In Theorem below we provide lower bounds for the error of any strong approximation of 
X^(T) for the processes X'^ from Section based on the whole path of (IT(t))tg[o,T] up to 
a time interval (to^^i) ^ [0 ,t'i/ 2]. The main tool for the proof of Theorem is the following 
simple symmetrization argument, which is a special case of the concept of radius of information 
used in information based complexity, see [29]. 

Lemma 1. Let P) he a probability space, let (f2i,^i) and (fl 2 ,-^ 2 ) be measurable spaces, 

and let Ih: —)■ fli and V 2 , I^', V!^-. —)■ 122 be random variables such that 

(V IP(Vi,V2) = ^{Vi,V') = PiViYa") • 

Then for all measurable mappings <h: f2i x 122 M o^nd (p: f2i —)■ M we have 

( 8 ) >lE[mv^,v')-^v^,v;%. 
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Proof. Observe that Q ensures that 

(9) E0<h(Ki,K2)-y^(V^i)|] 

This and the triangle inequality imply that 

(10) E[|<h(l/i,l/2)-¥^(lh)|] >lE[|<l>(lh,^')-$(l^i,14")l], 

which hnishes the proof. □ 

In addition, we employ in the proof of Theorem the following lower bound for the hrst 
absolute moment of the sine of a centered normally distributed random variable. 

Lemma 2. Let (O, A, P) be a probability space, let r G [1, cxd), and letY-.Q^Rbea r^)- 
distributed random variable. Then 


( 11 ) 

Proof. We have 


E[|sin(F)|] > 


■y/Svr 


■ exp( — 


TT 


( 12 ) 


E[|sin(y)|] £|sin(rz)|exp(^-y)d 2 ; 


> 


1 


This and the fact that 


TT 


expl—— ) / |sin(r 2 ;)| = 




expl — 


r 


TT 


TTT 

2 


|sin( 2 ;)| dz. 


2 


sin(a;)| dx > 


■W-f 


sin(a:)| dx = [rj 


sin(x) dx = [rJ > - 


□ 


(13) 

complete the proof. 

We hrst prove the announced lower error bound for strong approximation of X^{T) in the 
case of the time interval (toi^i) being sufficiently small. 

Lemma 3. Assume the setting in Section^ let ai, a 2 , a^, A, (3 G (0,cxd), and 7 G M &e given 
by 


r*ri 


(14) 


CKl = 


A = 


\f{s)\^ds, 02 = sup |/'(s)p, 03 = inf |/'(s)|^ 
se [0,0/2] se [0,0/2] 

rT 


mm 


ai 1 


2a2' (^2 


1/3 


r.T 3 


lfi'(s)l ds, 7 = / h{s)ds, 


'T3 


let xf G C°°(]R, (0, 00 )) be strictly increasing with hminf 3 ._,.oo V’()r) = 00 and = 1, let 

toAi G [0,ri/2] satisfy 0 < ti — to < A, and let u: (^([Ojfo] U —)■ M &e measurable. 


(15) E |At(T)-'t.((iy(i,)).,|„,,„,u|.,.T|) 


> 


87r3/2 


exp 


fj 


-1 


vT2 ^ 

(0-40)3/2757; 
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Proof. Define stochastic processes W, B : [to, ti] x D —)■ M and W ■. ([0, to] U [ti, Tj) x D —)■ M by 

(16) W(t} = ■ W(h) + . W(to), B(t) = W(t) - W(t) 

for t G [to,ti] and by W{t) = W{t) for t G [0,to] U [ti,T]. Hence, H is a Brownian bridge on 
[to,ti] and B and {W,W) are independent. 

Let Hi, 1^2: D —)■ M be random variables such that we have P-a.s. that 

( 17 ) 

ri= f° f{s)dW{s)+ r f(s)dW(s) + f(tt)W(t,) - f(to)WM - [' f'(s)W(s)ds, 

Jo J t\ J ^0 

Y 2 = - r f{s)B{s)ds 

Jto 

and put 

( 18 ) ai={E[\Yi\^]Y'^ 

for i G {1,2}. By the independence of B and {W,W) we have independence of Yi and l 2 - 
Moreover, for all t G (1, 2} we have Py- = 7V(0, af). Furthermore, Ito’s formula proves that we 
have P-a.s. that 


(19) Xt{n)=Yr + Y2. 

Therefore, we have P-a.s. that 

(20) Xt(T) = 7 . cos((ri + Y^)i,(Xt(T^))). 

First, we provide estimates on the variances |cri|^ and \o' 2 ^■ The fact that H is a Brownian 
bridge on [to,ti] shows that for all s,m G [to,ti] we have 

(ti — max{s, m}) ■ (min{s, u) — to) 


( 21 ) 


E[B{s)B{u)] = 


ih - to) 

In addition, the assumption infsg[ 0 Ti/ 2 ] |/Xs)| > 0 implies that for all s, m G [0,ri/2] we have 
/'(s) • /'(m) = |/'(s) • f'{u)\. The latter fact and (21) yield 

rii 


|(J2| = E 


( 22 ) 


f{s)B{s) ds 


l Jto 

'tl ftl 


'to Jto 


f{s) f{u)E[B{s)B{u)] dsdu 




(tl — max{s, u}) ■ (min{s, u} — to) 


Jto Jto 

Furthermore, it is easy to see that 

(tl — max{s, u}) ■ (min{s, u} — to) 

Jto J to 


{tl — to) 


ds du. 


(23) 


(<i - h) 


ds du = 


{tl - tp)-^ 
12 
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Combining (22) and (23|) proves that 

( 24 ) 


0(3 (ti — to)^ , I |2 ^ (^1 — toY 

U < -tt:- ^ *^2 S “ 


12 


12 


Next (24) and the assumption ti — to < A imply 
(25) 


CT 2 I < 0(2 |A| = min{Q(i/2,1} . 


By (19), by the fact that Yi and 1^2 are independent centered normal variables, and by (25) we 
get 

IcTi 

(26) 


|ai|'= E[|Fi| 2] =E[|Fi + F2p] -E[|F2p] -2 E[FiF2] 


= E[|X|’(ri)|^] - |(T2|^ = oi - \a2f > OLxj^l > |a2|^ , 
which jointly with (25) yields 

(27) |a 2 |^ < min{|(Ti|^ , 1 } . 

In the next step we put up the framework for an application of Lemma Observe that (20) 
and the assumption 7 7 ^ 0 imply 


(28) 


E 


\xt{T)-u{W)\ 


= hi -E 


|cos((yj + F2)^W(r3)))-^-u(LL)| 


Clearly, there exist measurable functions d**: C([0,to] U [fi,l],]R) —)■ M, i G {1,2}, such that 
we have P-a.s. that Yi = $i(hL) and X^{to) = $ 2 (hL). Moreover, by the independence of B 
and {W,W) we have independence of 12 and W. Therefore, we have = P^ ® Pyj = 

Pp^ ® P_y 2 = _Y^y We may thus apply Lemma with Oi = C([0,to] 0 [ti, 1],P), O 2 = M, 

Cl = W,V 2 = = ¥ 2 , VY = -I 2 , = ^ ■ M, and $: C'([0,to] U [fi,T],M) x M ^ M given by 

<I)(tn, y) = cos(($i(tn) + y) for w G C([0, to] 0 [fi, T], M), y e U to obtain 


E 


jcos((yi + U)iA(X*(r 3 ))) - i ■ u{W)\ 

= E 


(29) 


>i-E 

= i-E 


cos((<i)i(iy) + ¥ 2 ) y{^2{w))) - ip{w)\ 

cos((<I)i(W) + ¥ 2 ) Y{^ 2 {W))) - cos((<l>i(W) - ¥ 2 ) Y{^ 2 {W)))\ 
cos((Fi + F 2 ) Y{xt{rY)) - cos((yi - F 2 ) Y{Xt{A))) \ 

The latter estimate and the fact that Vx, ?/ G M: cos(a;) — cos{y) = 2sin(^^) sin(^^) imply 


E 


|cos((y + iyiA(x*(T3)))-y«(»')| 


> E 


|sin(yi?/>(x|'(r3))) ■ sin(y'2 ^’(^^(t's))) 


( 30 ) 
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Since {W, W), B, and (W(t) — W{T 2 ))t&[T 2 ,T 3 \ are independent we have independence of Yi, > 2 , 
and Xf{T 3 ) as well. Moreover, we have = 7V(0,/3). The latter two facts and (30) prove 


E 


jcos((yi + y) iix*{T3))) - i. «(»•) I 


(31) 

> [ E 

Jr 

sin('^(a;)Yi) 

■E 

sin (a;) 17 ) 


= [ E 

Jr 

1 sin('0(a;)yL) | 

■E 

1 sin('0(a;)F2) | 


P 




Next we note that (27) ensures that l/cr 2 > 1. This, the assumption that ip is continuous, the 
assumption that lim 3 ;_^oo ip{x) = oo, and the assumption that 1 p{^/^) = 1 show 


(32) 

It follows 


1/0-2 e [ip{\/^),oo) (ZipiR). 


(33) 


E 


> 


> 


sin('^(a;)Yi) 

f2tp-^{l/a2) 


■E 


sm 


{'ip{x)Y 2 )\ 7 ^ exp(-g) da; 


E 


'ip hi/f^2) 

1 


sin(-^(a;)Yi) 


■E 


sm{iP{x)Y 2 ) I ^ exp(-g) dx 


7273 






r-2i/) hl/'^2) 


E 


'ip 71 / 0 - 2 ) 


sm 


{il){x)Yi)\ 


■E 


sm 


{iIj{x)Y2)\ 


dx. 


We are now in a position to apply Lemma Observe that ( [2^ and the assumption that ip is 
strictly increasing imply that for all x G ['^“^( 1 / 0 - 2 ), 00 ), i G {1, 2} we have <Jiip{x) > aij 02 > 1- 
Employing Lemma we thus conclude that 

E[| sin(V;(a:)Yi)|] ■ E[| sin(V’(a;)y 2 )|] dx 


(34) 


'i/>“71/o-2) 


> 


-27 71/0-2) p -|2 1 „ 


'7-71/0-2) 


Furthermore, (24), (32), and the assumption that ip is strictly increasing ensure that 

(35) 


Vx/sy (ti-to)3/2;- 


Combining (31)-(35) proves 
E 


(36) 


cos((y + y)i/>t7(7-3))) - / ■ “(IT) 
1 


> 


\/27r/3 


exp 


(-||Y‘(^-^ir7|37j)| ) ■ 
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Finally, note that (32) and the assumption that "0 is strictly increasing imply 
Hence, we derive from ([M| that 


E 


(37) 


|cos((yi + Vi) i>{xt(n))) - i. u[W)\ 


>exp(-| 


-l( ^ 

a/^ (*1- 




□ 


This and (28) complete the proof of the lemma. 

We are ready to establish our main result. 

Theorem 4. Assume the setting in Section^ let Oi, 02 , Q^s, c, C G ( 0 , 00 ), and 7 G M &e 
given by 


(38) oi = / \f{s)f ds, 02 = sup |/'(s)p, 03 = inf |/'(s)p, /3 = j \g{s)f ds, 

Jo se[0,n/2] s6[0,n/2] 

rT 


•T2 


(39) 


7 = 


h{s) ds, c = 


I 7 I 


8 7r3/2exp(^)’ 


C = 


a/I2 max{l, T3/2y7f2} 


^a3min{l,y^} 

let -0 G C“(M, (0, 00 )) be strictly increasing with liminf 3 ,^ 00 '^(a^) = 00 and '0(\/2^) = 1, 
let 0 < to < ti < T 1 I 2 , and let u\ C([0,fo] U [ti,T],M) —)■ M &e measurable. Then [C'/(fi — 
to)^^‘^,oo) C and 

X^(T) -M((W(s)),6[o,to]u[ti,T])| > c-exp(^-| ■ 


(40) 


E 


Proof. Let A G (0, cxd) be given by (14). 

First, assume ti — to < A. By Lemma and by the properties of if we then have 


(41) 
and 

(42) 


c'/’(») 


E 


Xt{T)-u{{W{s)),^^oMAt.p)\ 


> c ■ exp 


if 


-1 


v42 ^ 




It remains to observe that 
(43) 


^/Vl 


< 


c 


and that if~^ is strictly increasing to obtain the desired result in this case. 

Next, assume that fi — to > Then Lemma together with the properties of if yield 

(44) [^3^^,cx)) C V^(M) 

and 


(45) 


E 


|.Y*(r) - «((w^(s)).5[„,„iu[,.,T]) 


> c■exp 


A 


if 


-1 / F 12 

I A3/27«J 
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Since 

(46) 






< 




J^3I2 


< 


c 


- v^min{l, {h - tof/^ " (t, - to)"/' 

and since is strictly increasing, we obtain the claimed result in the actual case as well. □ 

Theorem [^implies uniform lower bounds for the error of strong approximations of the solution 
processes X'^ in Section at time T based on a finite number of function values of the driving 
Brownian motion W. This is, in particular, the subject of the following corollary. 

Corollary 1. Assume the setting in Section^ let Oi, 02 , cis,/S, c, C G (0, 00 ), and 7 G M &e 
given by 


(47) 

(48) 


bTi 


(y.\ — 


7 = 


|/(s)|'ds, 02 = sup |/'(s)|', 03 = inf |/'(s)p, (3 = 

se[0,n/2] s6[0,n/2] 

rT 

h{s)ds, c= --i-t-'- C = 




\g{s)f ds, 

■\/l2 max{l, T^/^yTf2} 


'T2 


'ra 


8 7r3/2exp(^)’ 


4^ A/a3min{l, 

and let ip G C°°(M, (0, cxd)) be strictly increasing with liminf 3 ,^ 00 '^(t) = 00 and '0(\/^) = 
1. Then for all n G N n [2T/ri,oo) and all measurable u: C{[T/n,T],'R) —)■ M we have 
[Cn^/^T”^/^, cx)) C i/’(M) and 

(49) E \xf{T) -u{{W{s))s(.[T/^,T])\ > c-exp[-| ■ -n"/')'^ 


for all n eN, Si,..., Sn E [0, T] and all measurable u: 
-(/^(M) and 


-E 


we have [ 8 C'n^/^(ri) ^/^, 00 ) c 


(50) 


E 


Xf{T) M(iy(si),. . . , hh(Sn)) ij > c • exp 


and for all n E Nn[2T/ri, 00 ), si,..., G [0, T] and all measurable u: M”xC([T/? 7 ,, T], M) —)■ M 
we have [2^/^ C ■ n"/T^/^, 00 ) C and 

(51) E [|X 4 ^(T) - u{W{s ,),..., W{sn). (PB(s)),e[T/„,T]) |] > c ■ exp(-| • ■ n^) T) • 

Proof. Let n G N with T/n < ri/2 and let u\ C([T/? 7 ,,T],M) —)■ M be a measurable mapping. 
Then Theorem 0 with to = 0 and ti = T/n implies [C ■ n^/^/T^/^, 00 ) c ^jJ(R) and 

(52) E[|W*(T)-f(ar(s)).,|x/.,Tl)|] >c.exp(-|.p-'(37§57i)r)- 
This establishes (49). 

Next let n G N, si,... ,s„ G [0,T] and let u: —)■ M be a measurable mapping. Then 

there exist sq, si,..., s„+i G [0, T] such that 0 = sq < si < ■ • • < s^+i and {sq, si,..., Sn+i} 3 
{si, ... ,Sn, ti/ 2}. In particular, there exists i E {1, 2,..., n + 1} such that 


(53) 


Si < f and Si - Si_i > 
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Using Theorem 0 with to = Sj_i and U = Si and the fact that ^ is increasing we conclude 
that oo) c [C( 2 (n + oo) c [(^/(si — cxd) C and 


E 


(54) 


\Xt{T) - u{W{So), W{si), W{Sn), W{Sn+l)) \ 

>c-exp(-|- >c-exp(-|- 


This implies (50|). 
The proof of 


is analogous to the proofs of (49) and (50). 


□ 


In Lemma 1^ below we characterize a non-polynomial decay of the lower bounds in (49), (50), 
and (51) in Corollary in terms of a exponential growth property of the function 'ifj. To do so, 
we recall the following elementary fact. 

Lemma 4. Let <^i : M — )• [0, cx)) be non-decreasing, let ip2'- ^ ^ [0, C)o) be non-increasing, and 
assume that liminff^g^^oo + 1)] = oo. Then liminfjRga-^oo [</5i(3^) • T 2 {.x)] = oo. 

Proof. By the properties of ipi and ip 2 we have for all x G M that (pi{x) ■ (p 2 {x) > (^i([xj) ■ 
</ 92 ( [xj -|- 1). Hence 

(55) 


liminf [</ 3 i(x) • <^ 2 {x)\ > liminf [(pi{n) ■ (p 2 {n 1 )] = oo, 

NBn^oo 


which completes the proof. 


□ 


Remark 3. We note that in general it is not possible to replace in Lemma the assumption 
liminf \g)i{n) -(^ 2(11 + 1)1 = 00 by the weaker assumption liminf YPiiji) ■ g) 2 ip)\ = 00 . Indeed, 

NSn—>-00 ^ ^ N 9 n—>-00 

using suitable mollifiers one can construct (^ 1 , <^2 ^ C°°(]R, [ 0 , cx))) such that (^1 is non-decreasing 
with VnGZVxG [? 7 ,,n-f- 1/2]: (pi{x) = exp((n -|- 1/2)^) and such that ip 2 is non-increasing 
with 'i n ^TL'i X ^ [n — 1/2, n]: (^ 2 ( 2 ^) = exp(—n^). Then 


(56) 


liminf [.piiji) ■ g:>2ip)] = liminf exp((n -|- 1 / 2 )^ — n^\ = 00 , 

NBn^oo NBn^oo ^ 

liminf Vpi{n) ■ g:>2{n -\- 1 )] = liminf exp((n -|- 1 / 2 )^ — {n -\- 1 )^) = 0 , 

N9n^cxo NBn^oo ^ 

liminf [(pi{x) ■ ^ 2 {x)\ < liminf [(pi{n + 1 / 2 ) • (p 2 {n + 1 / 2 )] = 0. 

R9a:^oo NBn^oo 


Lemma 5. Let rii,7]2,r]3 G (0, 00 ) and let M —)■ (0, cx) be strictly increasing and continuous 
withWm.m.ix^aof^ix) = 00 . ThenM q G (0, x): liminf^g^^oo [n'^-exp(— 771 )] = 00 

if and only if\/qE (0, x) : liminf]R 3 a,^oo [t/(x) • exp(—gx^)] = x. 

Proof. We use Lemma with c^i(x) = x'^ and (f 2 {x) = exp(— 771 \'4)~^{'q2X^^)f) for x G M to 
obtain 

Vg G (0, x): liminf ■ exp( — 7/1 {ri2n^^)\^ )] = x) 

N9n—>-oo ^ / 

fVg G (0, x): liminf Tx'^ • exp( — 7 ]i {ri2x'^^)'\^ )] = xV 


(57) 
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Furthermore, 


Vge(0,cxD): liminf ■ exp( — ? 7 i 1-0 ^)] = oo 
fVg e (0, oo): liminf • exp( — r]i {t]2x'^^)\‘^ )1 = oo 

( 58 ) (Vg e (0, cxd) : liminf Tx'^ ■ exp ( —r/i |'0“^(772 t) I ^ )1 = oo) 

V RSa;—>-00 ^ ' I / J J 

fVg e (0, oo): liminf Tx • exp( — ^ |'0“^(r72x) 1^ )1 = oo) 

^ (y q G (0, oo): ■ exp( — ^ |'0“^(x)|^ )] = ooj. 

Using the properties of "0 we have 

fVg G ( 0 , oo): liminf Tx ■ exp( — ^ {x)\‘^ )1 = oo) 

( 59 ) ^Vg G ( 0 , 00): hminf [■^(x) • exp( — ^x^)] = ooj 

■v^ (Vg G ( 0 , 00): liminf ['^(x) ■ exp( — gx^)l = 00), 

which completes the proof. 


□ 


As a immediate consequence of (51) in Corollary and Lemma we get a non-polynomial 
decay of the error of any strong approximation of X^{T) based on n G N evaluations of the 
driving Brownian motion W and the path of W starting from time T/n if -0 satisfies the 
exponential growth condition stated in Lemma 

Corollary 2. Assume the setting in Fecfzonj^ let (3 G (0,cx)) he given by [3 = fr' \ 9 {s)f ds, 
and assume that 0 G ^“(M, (0, 00 )) is strictly increasing with the property that = 1 

and\/q G (0,oo): liminfa,_,.oo [ipix) • exp(—gx^)] = 00 . Then for all q G (0, 00 ) we have 

(60) 


lim inf ( ■ inf 

n^oo y u: R‘^xC{[T/n,T],R)^R 
measurable, 5i,...,Sn£[0,T] 


E 


\xt{T) - n(lU(si),..., Wisn), (lU(s)),e[T/„,r]) 


= CX). 


The following result shows that the smallest possible error for strong approximation of X^(T) 
based on n G N evaluations of the driving Brownian motion W and the path of W starting 
from time T/n may decay arbitrarily slow. 

Corollary 3. Assume the setting in Fecfionj^ let (3 G (0, 00 ) he given by (3 = |g(-s)|^ ds, and 

let (a„)„gN C (0, 00 ) satisfy hmsup„_,.f^ a„ = 0. Then there exist a real number k G (0, 00 ) and 
a strictly increasing function 0 G ^“(M, (0, cxd)) with liminfa-^oo 0 (t) = 00 and 0(a/^) = 1 
such that for all n G N, Si, S 2 , ..., Sn ^ [0, T] and all measurable u: M” x C[[T/n,T],M/j —)■ M 
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we have 
(61) 


E 


\Xt{T) - u{W{s^), . . . , WiSn), {W{s))s^^T/n,T]) \ 


> K - Ur,. 


Proof. Without loss of generality we may assume that the sequence is strictly decreas¬ 
ing. Let c, C e (0,oo) be given by (39) and put C = Choose no G N fl P^/ti,cx)) 

such that for all n G {no, no -|- 1 ,... } we have 

( 62 ) On < I < C ■ and | > 2/?; 

and let (&n)ne{no-i,no,...} C ( 0 , oo) be such that 6^0-1 = cind such that for all n G (no, no -|- 
1,... } we have 

.1 1/2 


(63) 


^77, - 


IMP 


Note that {bn)n£{no-i,no,...} IS strictly increasing and satishes hm„_,.oo bn = 00 . 

Next let -0: M —)■ (0, cxd) be the function with the property that for all n G (no, no -l- 1,... }, 
a: G M we have 

1 , 


(64) fj^x) = < 


1 + 


C ■ nQ — 1 


C-n^, 


C-{n-lf + 


C-n^-C-{n-lf 


if X < bno-l, 
if X = bno-l, 

if X G {bno-l,bno), 

if X = bn and n > uq, 

, if X G {bn-i, bn) and n > uq. 


Then is strictly increasing, positive, and inhnitely often differentiable and V’ satishes V’(\/^) = 
1, liminfa;_,.oo '0(2^) = oO) and = (0, cxd). 

In the next step let Sn G [0, 00 ), n G N, be the real numbers with the property that for all 
n G N we have 


(65) En = inf inf 

si,...,s^G[ 0 ,T] u: M^xC([T/n,T],]E 
measurable 


E 


1X4^(T) - n(W(si),..., W{sn). (W(s)),er/„,r])| 


Estimate (51) in Corollary yields that for all n G {no, no -|- 1,... } we have 

(66) En > c-exp(^-| ■ = c-exp(^-| • 

Since the sequence (£n)neN is non-increasing, we have for every n G {1,2,... ,no} that En > 
Eno > c - a„g. We therefore conclude that for all n G N we have 


C • dyi' 


(67) 


ca 


En>c- min{l,a„(,/a„} ■ ^ ■ a„. 


Oi 


which completes the proof of the corollary with k = c ■ an^/ai. 


□ 
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Next we extend the result in Corollary to approximations that may use hnitely many 
evaluations of the Brownian path as well as the whole Brownian path starting from some 
arbitrarily small positive time. 

Corollary 4. Assume the setting in Section let (3 G (0,cxd) he given by 13 = fr2 
and let (a„)„gN C (0,cx)) and {6n)nen C (0, cx)) satisfy lim„_^oo On = hm^^oo = 0. Then 
there exist a real number n E (0, oo) and a strictly increasing function -0 G C°°(]R, (0, oo)) with 
hminf 3 ;^oo fj{x) = oo and = 1 such that for all n E N, si, S2, ..., G [0, T] and all 

measurable u: MT x C([5„,T],M) —)■ M we have 


( 68 ) 


E 


\X*{T) - '(.(IT)*!),.... W{s„), (ir (s)),e|j„,T|) 


> K ■ Or 


Proof. Without loss of generality we may assume that the sequence (5„)neN is strictly decreasing. 
Let (fcn)neN C (0, oo) be the strictly increasing sequence of positive integers with the property 
that for all n G N we have 


(69) kn = FAn] + n. 

Moreover, let (a„)„gN C (0, oo) be a sequence such that for all n G N we have 

( 70 ) On 

and lirnm^oodm = 0. Then Corollary implies that there exist a real number k E ( 0 ,oo) and 
a strictly increasing function G C“(M, (0,oo)) with liminf 3 ._^oo A(2^) = oo and '0(\/^) = 1 
such that for all n G Si,S 2 ,... ,Sn G [0,T] and all measurable u: M"' x C(AA,T],M) —)■ M 
we have 


(71) 


E 


\xt{T) - u{W{s,), ..., W(sA, (lT(s)).e[T/n.r]) | 


> K - a„. 


Let n G N, let m: X C'([5n,T],M) —)■ M be a measurable mapping, and let si,S 2 , .. .,Sn ^ 

[0,T]. Note that (69) implies 6n > '^/k„ and kn > n. Put Sm = for m G {n + l,n + 

2,..., kn}. Clearly, there exists a measurable mapping u : x C([T/fc„, T], M) —)■ M such that 

w(W(si),...,W(sA,(W^(A).e[^„,T]) = w(lT(si),...,W(sfcJ,(W(s)).6 [T//c„,T]). Hence, by (0 
and by ([70|), we have 


(72) 


E 


\xtiT) - qw')®.),.... W{s„). (M 4 s)).,[j.,ti) I 


> K ■ ak„ = K ■ On, 


which completes the proof. 


□ 


5. Upper error bounds for the Euler-Maruyama scheme 

A classical method for strong approximation of SDEs is provided by the Euler-Maruyama 
scheme. In Theorem below we establish upper bounds for the root mean square errors of 
Euler-Maruyama approximations of X^{T) for the processes X^, ip E C'°°(M, (0, cx))), from 
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Section In particular, it turns out that in the case of non-polynomial convergence the Euler- 
Maruyama approximation may still perform asymptotically optimal, at least on a logarithmic 
scale, see Example below for details. 

We first provide some elementary bounds for tail probabilities of normally distributed random 
variables. 


Lemma 6. Let (hi. A, P) be a probability space, let x G 
normal random variable. Then 


and let Z: Q ^ M. be a standard 


(73) P(Z>x) < ^•exp(-^). 

Proof. For every i/ G [0, cx)) we have 

(74) {y + xY - x\x\ - Y = “ I^D) = Uy'^ + '^^y + 4a:^l(-oo,0](a;)) > 0. 


Hence 


(75) 


nz>x)=l ^■exp( 


jy+P 


j dy 


poo 

<exp(-*l) f ^.exp(-^)dy=i.exp(-^). 


which completes the proof. 


□ 


Lemma 7. Let (hl,^, P) be a probability space, let a G [0,cxd), c G (0, cxd) fl [a, cxd), and let 
Z: Q ^ M. be a A/'(0, a"^)-distributed random variable. Then for all x eM. we have 


(76) 


P(Z > x) < exp( - 




Proof. In the case a = 0 we note that for all a; G M we have 
(77) P(Z >x)= l(_„,„|(i) < exp ( - 


2c2 


)■ 


In the case a > 0 we use Lemma to obtain that for all x G [0, oo) we have 
(78) P(Z >x)= P(f > p < i .exp( - ^) < exp( - fl), 

which completes the proof. 


□ 


Next we relate exponential growth of a continuously differentiable function to exponential 
growth of its derivative. 

Lemma 8. Let G C^(M,M) satisfy'iq G (0,oo): hminfa;p^.oo ["^(a;) • exp(—ga:^)] = oo and 
assume that xjj' is non-decreasing. Then Vg G M: liminfa^p^oo Wi.x) ■ exp(—ga;^)] = oo. 

Proof. Since Vg G (0, oo): liminf 2 ._ 5.00 [f’ix) ■ exp(—ga;^)] = 00 , we have 

(79) VgGM: liminf ['^( x) ■ exp(—gx^)] = 00 . 
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By the fundamental theorem of calculus and the assumption that rjj' is increasing we obtain for 
all X G (0, oo) that 


(80) 


^ 1 r ^, 1 r ^ 

ijj{x) = - ilj{x)dy>- 7p{y)dy=- 


X 


X 


Hence, for all g G M we have 

liminf [^jJ'{x) ■ exp(—> liminf 


(81) 


'0(x) — -0(0) 

X ■ exp(ga;2) 


X 


> lim inf 

X^OO 


^jj{x) - l^/J{x) 

X ■ exp(ga;2) 


= liminf 


'ip{x) 


x^oo \_2x ■ exp(ga;2) 
which completes the proof. 


> liminf Ujj{x) ■ exp(—2ga;^)l = oo, 


□ 


We turn to the analysis of the Euler-Maruyama scheme for strong approximation of SDEs in 
the setting of Section 

Theorem 5. Assume the setting in Section^ assume that Ti < T 2 , let [3 G (0,oo) be given 
by ^ = Ir' \g{s)f ds, let S G (0,1), let ^|J G C°°(M, (0, oo)) he strictly increasing such that 
= 1, such that Vg G (0, oo): hminfa,,_^.oo ["0(2:) • exp(—gx^)] = oo, and such that tjj' 
is strictly inreasing, and let : {0,1,..., n} x —)■ IR"^, n G N, satisfy for all n G N, 

fc G (0,1,..., n — 1} that = 0 and 


(82) = iff-"’ + ^♦(if“•">) I + <t(V»'”') 

Then there exist real numbers c G (0, oo) and no G N such that [|r 
that for every n G (uq, Uq + 1,... } we have 


io|'^,oo) C fj'iR) and such 


(83) fE[||X^(T) 


l|2 

\ 1/2 r 

IIr4 

) 


exp 




+ exp 


/\-i/ 


n 


')f 


Proof. Throughout this proof let APE”: fl —)■ M, j G (1,2,... ,n}, n G N, be the mappings 
with the property that for all n G N, j G (1, 2,..., n} we have AIE” = W{^) — let 

/!„ G M, n G N, and 7 „ G M, n G N, be the real numbers with the property that for all n G N 
we have 


(84) 


'ln = J2n-K 

i=i 


(j-i)r 


). 


(3n ^ 

1=1 


9{ 


jj-TT 

n 



and let : (0,1,..., n} x fl —)■ M, / G (1,2, 3,4}, n G N, be the stochastic processes with 

the property that for all n G N, fc G (0,1,..., n} we have ..., x|’^’"'^). By the 

properties of f,g, h stated in Section]^ and by the dehnition of yA and a (see (|^), we have for 
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all n G N, A; G { 0 , 1 ,..., n} that ^ and 




mm{k,\nri/T ~\} 


k 


•hr’=EsirF) ■ 


i=i 

min{ k, \nr^ /T]} 


J-! 


(85) 


i=i 

k 


j=[nT2/T\+2 


X 


{'ip,n) 

4 ,fc 


E r rrr) • ■ v'(Airi)) 


i=i 


= E rM“G^)-<=o=(rri'V'(Ar-”i)). 

j=[nT3/rj+2 

In particular, for all n G N, fc G [^,cxd) fl {1,2,..., n} we have and for all 

n G N, /c G [^, cxd) n {1, 2,..., n} we have = X^^\ Therefore, for all n G N we have 

(86) Air’= E pMrr)'“hAir’-v4Air’))=7n'cos(.?r”’'<f(.hr’))- 

j=[nT3/TJ+2 

We separately analyze the componentwise mean square errors 

(87) T,„ = E[|A*(T)-rr’l'] 

for i G {!,... ,4}, n G N. Clearly, for all n G N we have Si^n = 0. Moreover, Ito’s isometry 
shows that for all n G N we have 

( 88 ) 


^ 2 ,n — IE 


f>jT/n 


E/. (/(i>) -/(“lE)) <'‘r(i>) 


j = l J{j-l)T/n 


riT/n 


E/, i/(*)-/(rE)r* 


J(j-l)T/n 


n « 

< sup i/'(i)rE / 

lefo.nl Jo 


ie[o,Ti] 

and, similarly. 


'■jT’/n ry-i3 

te[o,ri] 


J = 1 J(i-l)T/n 


(89) 


r3 


£ 3 ,n < ^ ■ sup |^'(^)|^ < T- sup |/(t)|^. 


3n2 


le[T2,r3] 


te[0,ri] 


(14)). From ( 86 ) we obtain 
|2 


We turn to the analysis of £ 4 ^^, n G N. For this let 7 G M be given by 7 = h{s) ds (see 

+ 2|7-7n|^ 


(90) £4,n < 2 |7|^ ■ E I cos(W 2 ^(T) • ^{Xt{T))) - cos(x(P’") • 


-(■?/',[2 
l,n 
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Clearly, for all n G N we have 


(91) 



^ ■ sup |ft'(()|. 

ie[T3,T] 


Using a trigonometric identity, the fact that Vx G M: |sin(x)| < min{l, |x|}, inequality (88), 
the fact that = A/'(0, /3), a standard estimate of Gaussian tail probabilities, see, e.g., [T71 

Lemma 22.2], and the fact that = VW we get for all n G N that 


E 


(92) 


j cos(X*(r) . i,(X*(T))) - cos(.?'V> . i'(X*(T))) 

= 4 . E I sm(l (X*(T) - .?<*•"') iiX*(T))) siu(l (X*{T) + x‘'*f) V-(V*(T))) 
< 4 . E I sm(l (X*(T) - <l,(X*(T)) 

X*{T) - |,i.(At(T)) f 


< E 

< n^'^E 


< 


3n2(i- 


+ 4-P(^X3^(T) > 

|.V(T) - At>r] + ■ exp(-A . 


By Lemma we have Imvx^oo {x) = oo. Hence, there exists ni G N such that [|ni|'^,oo) c 
'0'([O, cx))). Put 


(93) 


no = max< rii 




and let n G {no, no + 1,... }. Then {W{s))s(.[Q,\nT-,/T\T/n] and (lU(s) - lU(r 2 ))se[r 2 ,r] are inde¬ 


pendent, which implies independence of the random variables and Xj(T) Using 

the latter fact as well as the fact that is strictly increasing and the estimates in (89), we may 


:(p,n) 
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proceed analoguously to the derivation of 


E 


to obtain 

('?/’,n) / / 0(n)\ \ |2 


(94) 


|cos(.?'*'”' . <f(.Y*(T))) - cos(A'g’”> . V'(.?S))| 

< 4 ■ E[|sm(l. A»„”> . [i-{X*{T)) - V'(A»„”>)]) 

< E [b£"> r ■ |V-(A*(r)) - 

+ 4.P(4''(max{A3*(T),A“„">}) >n') 


< ■ E 


X, 


(' 0 ,n) 12 
2,n I 


■E 

-b 


|.Y 3 *(t) - X 


12 


3,n 


< 


+ 4.P(max{X3’’(r),X«;"'} > (V/)-‘(ii")) 

■ sup |/(i)r sup |g'(«)p 

tS[0,ri] lS[T 2 ,r 3 ] 


+ 4.P(A'3*(r) > (V/)-‘(n‘)) +4-P(X<*-”> > (V/)-‘(u‘)). 

Note that = A^(0,/S) and = A/'(0,/dn) and sup^gj^/3m G [/3, oo). We may therefore 

3 3,n 

apply Lemma to conclude 

P(X3*(T) > (V.')-‘(n‘)) +P(xg'”' > W)-\A) 

|(P')-l(nb|^ , l(h')-HnbP' 


(95) 


< exp(- 


-) + exp(- 




Combining (90)-(92) and (94)- 


2/3 / nv 2sup„3gf,/3„i . 

ensures that there exist ci,C 2 G ( 0 , oo) such that for all 


n G {rio, rio + 1 ,... } we have 

(96) £4 < Cl ■ + exp(-C 2 • +exp(-C 2 • |(i/’')"^(^'^)T)) • 

By assumption we have for all q G (0,oo) that liminf^i^oo \^{x) ■ exp(—ga;^)] = oo. Hence, 
Lemma ensures that there exists C 3 G (0, 00 ) such that for all n G N we have 


(97) 


( 1 ^ < C 3 ■exp(-C 2 1-0 ^(n^)f) 


Combining ( 88 ), (89), (96), and (97) hnishes the proof. 


□ 


Example 2. Assume the setting in Section]^ assume that ri < r 2 , let /3 G (0, 00 ) be given by 
/3 = |g(s)|^ds, let '0;: M —)■ (0,oo), I G {1,2}, be the functions such that for all a; G M we 

have 


(98) 

(99) 


= exp(x3 + 2 a: - - 2{2l3f/^) , 

' 4 ) 2 {x) = exp(a:exp(x^ + l) “ (2/3)^/^exp(2/3 + 1 )) , 
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and for every n G N, / G {1, 2} let : {0,1,... , n} X —)• be the mapping such that 

for all fc G {0, 1, 2, ..., n — 1} we have = 0 and 

(100) ^ + a(xf'’"^) (hh(^^) - fh(f)). 

Clearly, we have '0i)V’2 £ (0, cx))) and ■^ 1 (^ 2 ^) = 'ip 2 {\/W) = 1- Moreover, for all 

q G (0, cxd) we have 

(101) liminf \'iIji{x) ■ exp{—qx‘^)] = liminf [' 02 ( 2 ^) ■ exp(—gx^)] = 00 . 

a:i-^oo ^ a::i-^oo ^ 


Furthermore, for all a: G M we have 


( 102 ) 


'0i(a:) = (3x^ + 2) • 'il^i{x) > 0, 

^jJl{x) = [6x + (3a:^ + 2)^) • ^jJl{x) = (Ox^ + 3a;^ + 3 + (3a: + 1)^) • ^pl{x) > 0 


and 

= (2a:^ + 1) exp(a:^ + l) • ‘ip 2 {x) > 0, 

'ip 2 {x) = (4a: + (1 + 2a:^) (2x + (1 + 2x^) exp(x^ + l))) exp(x^ + l) ip 2 {x) 

^ ^ > (4a: + (1 + 2a:^) (2a: + 2(1 + 2a:^))) exp(x^ + 1 ) t/’ 2 (a:) 

> (4a: + 7/4 ■ (1 + 2a:^)) exp(a:^ + l) 'ip 2 {x) > (17/28) exp(x^ + l) 'ip 2 {x) > 0. 


Hence, 'ipi, 'ip[, 'ip 2 , and ■02 strictly increasing and we have = V’ 2 (^^) = (O 5 '^)- 

Using Corollary and Theorem with 5 = 1/2 we conclude that there exist ci,C 2 G (0, cx)), 
no G N such that for all k G {1, 2} and all n G {no, no + 1, • • • } we have 


(104) 


ci-exp(-|- 


< (e 


y/2 

VI 



Next, we provide suitable minorants and majorants for the functions k G (1, 2}, and 

(V’{)~^ k G {1,2}. To this end we use the fact that for all a G M and all strictly increasing 
continuous functions /i, / 2 : [a, 00 ) —)■ M with fi > f 2 and limj-^oo f 2 {x) = 00 we have 

Mx € [/,(a),oo): x = Mf^'ix)) < fidi^x)) 


(105) 

and therefore 

(106) 


Vx G [/i(a),cx): //(t) </s ^(a:). 
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Clearly, for all x G [1, cxd) we have 

(107) 


exp {x^ + 2 — (2/5)3/2 _ 2(2/?)1/2) < 

< exp(3a:^) , 

exp(exp(x^) — (2/3)^/^ exp(2/3 + 1)) < ' 02 (a;) < exp(exp(a;^ + a; + l)) < exp (exp (Sx^)) , 
'ip'i{x) < exp(3a:^ + 2) • '4)i{x) < exp(8a:^) , 

V> 2 (ai) < exp(3a:^ + 2) • 'il) 2 {x) < exp (exp (Sa:^)) . 


We may therefore apply ( 106| ) with a = 1 to obtain that for all x G [exp(exp( 8 )), oo) we have 

(ln(a:) - 2 + {2/3f^^ + 2(2/3)^/^)^^^ > > 3"^/^ • (ln(a;))^/^ 

(ln(ln(a;) + {2/3^^^ exp{2/3 + 1)))'/' > ij^\x) > 3-'/' ■ (ln(ln(a;)))'/2 

> 8-^/3 . (ln(a;))V3, 

(V'2)~^(^) ^ 8“^/^ ■ (ln(ln(x)))^/^. 


(108) 


Combining (104) with (108) shows that there exist Ci, 02 , 03,04 G (0, 00 ), rio G N snch that 


for all n G {no, no + 1,... } we have 

1 /2 

0i-exp(-02-|ln(n)p/3) < (^e[||X^H2^) - ) < 03 ■ exp(-04 ■ | ln(n)p/3), 


(109) 


Oi • exp (-02 ■ ln(ln(n))) < (e\\\X^^{T) - 


1/2 


< 03 • exp (—04 ■ ln(ln(n))). 


In particnlar, in both cases the Enler-Maruyama scheme performs asymptotically optimal on a 
logarithmic scale. 


6. Numerical experiments 

We illnstrate our theoretical hndings by numerical simulations of the mean error performance 
of the Euler scheme, the tamed Euler scheme, and the stopped tamed Euler scheme for a 
equation, which allows a decay of error not faster than 0 ■ exp( — i/c ■ | ln(n)p/^) in terms of 
the number n G N of observations of the driving Brownian motion, where c G (0, cx)) is a real 
number which does not depend on n G N. 

Assume the setting in Section]^ assume that T = 1, ti = r 2 = V 4 , T 3 = 3 / 4 , assume that for 
all a; G M we have 

fix) = l(-oo,i/4)(a:) • exp(^31n(10) + 
g{x) = 1 (i/ 4 , 3 / 4 )(x) ■exp(^ln(2) +41n(10) + - ^ 3 ^ 74)1 

h{x) = 1 ( 3 / 4 , 00 ) (ai) • exp(^41n(10) + 37 ^), 


(110) 
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(cf. Example 1), let (3 G (0, oo) be given hy (3 = \g{s)f ds, and let V’: t (0, oo) be the 

function such that for all x G M we have 

^jJ{x) = exp(x^). 

Recall that the functions /, g, h, and 'ip determine a drift coefficient —)■ and a 

diffusion coefficient cr: —>■ see Furthemore, recall that the fourth component of the 

solution X'3 of the associated SDE at time 1 satishes that it holds P-a.s. that 


( 111 ) 


Xt{l)= f h(s)ds-cos(Tf(s)dW(s)-P;lf^^;^g(s)dlV(s) 

3/4 Vo ^ ' 


f3/4 


see 


(§. 

Furthermore, let Aj//Ai//Ai”’/) : {0,1,..., n} X fi ^ R^ ii £ N. 

g G {1, 2, 3}, be the mappings such that for all g E {1,2, 3}, n G N, k E {0,1,..., n — 1} we 
have = 0 and 

-?h/ = At’-' + yXX) i + £(.?h’‘) (‘r( A) - ‘r(l)), 


j^{n ),2 _ j^{n ),2 ^ 




- + oAA')(»r(^)-iF(p), 


( 112 ) 


1 + i 


A(ii),3 _ ^(n),3 

-^k+1 ~ 


+ 1 


{lXfe"'’^llR4<exp(|ln(n)|V2)| 


^0(jf(«).3) i ^ o(At>4) (W.-(/tl) - W-’d)) 


.1+iih*(AA^") /+(W'(¥) - W'’(f) 

Thus XW’i, X (")’3 are the Euler scheme (see Maruyama ng), the tamed Euler scheme 

in Hutzenthaler et al. [H], and the stopped tamed Euler scheme in Hutzenthaler et ah [TB] . 
respectively, each with time-step size 1/n. 

Let E [0,oo), n E N, g E {1,2,3}, be the real numbers with the property that for all 
n E N, g E {1, 2, 3} we have 

£j = E[|.Y*(l)-.?<Al]. 

let /: M — )■ M and -0: M —>■ ( 0 ,cxd) be the functions such that for all x G M we have /(x) = 
exp((2/3)^/^) • /(x) and ■^(x) = exp(—(2/3)^/^) • //’(x), and let cxi, 02 , 03 , c,C E (0, cx)) be the real 
numbers given by 

(113) tti = / |/(s)pcis, 02 = sup |/'(s)p, 03 = inf |/'(s)p, 

0 se[0,n/2] se[0,n/2] 


(114) 


c = 


lX- 3 ^(g) ds\ 

8 7r3/2exp(^)’ 


C = 


a/12 max{l, 
ya^min{l, ^/f} 


In the next step we note that ip E C°° (M, (0, 00 )) is strictly increasing, we note that lim infa;_,.oo 'V’(x) 
00 , and we note that 'tp{-\/2(3) = 1. We can thus apply inequality (50) in Corollary (with the 
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functions /, g, h, and 'ip) to obtain that for all n G N, Si,...,Sn G [ 0 , 1 ] and all measurable 
n: —)■ M we have cx)) C V’(M) and 


(115) 


E 


Xt{l)-u{W{si),...,W{sn))\ > c-expf-| • |^/> 


This and the fact that Vy G ■^(M): -ip ^{y) = [ln(|/■ exp((2/3)^/^))] ensure that for all n G N, 
si,..., Sn G [0,1] and all measurable n : M” — )■ M we have 


( 116 ) 


E 


At(l) 


>f.exp(-|.|ln(j^.exp((2Dp).n’/^) 



= f. exp(-| |in( «g“Piyi ) + I lri(n)p''’) 
> c. exp(-| |ln( »-^7g,y V ) 


ln(n)P^V 


In particular, this proves that there exists a real number c G (0, cx) such that for all g G {1,2, 3}, 
n G M we have 


(117) el = E[|At(l) - > c. exp( - P I ln(n)|"'’). 

In the next step let m = 5000, N = 2^^, let B = {Bi,... ,Bm)'- [0,1] x hi —)■ be an 

m-dimensional standard Brownian motion, and let = {Y ^,..., 1^): hi —)■ M, iV G N, be 
the random variables with the property that for all iV G N, /c G {1, 2,... ,m} we have 


(118) 




/ h{s) ds ■ cos 

^ 3/4 



Liv/4J 


E ni) ■ Mi) ■ i’ 


N 


L3iv/4J 

E 


i=\N/ 4 :'] 



The random variables Yj^, k G {1,2,... ,m}, iV G N, are used to get reference estimates of 
realizations of X^l)- Onr nnmerical resnlts are based on a simulation 


(119) (^li • • • 1 bm) — ((&l,i)ie{0,l,...,w}, • • • , {bm,i)i&{0,l,...,N}) £ 

of a realization of ((5i(yAr))jg{o,i,...,Ar},..., (5m(yAf))ie{o,i,...,v}) (a realization of (5i,..., Bm) 
evalnated at the eqnidistant times ^/n, i G (0,1,..., N}). Based on (6 i,..., bm) we compnte 
a simulation (yi,..., ym) G MB’ of a realization of (Y^ ,..., Y^) and based on (6 i,..., bm) we 
compnte for every 77 G (1, 2, 3} and every n G (2°, 2\ ..., 2^®} a simulation ..., G 

of a corresponding realization of m independent copies of Then for every g G (1, 2, 3} 

and every n G (2'^, 2^,..., 2^®} the real nnmber 


( 120 ) 


-| /ft 

m 

(.=1 


X 




serves as an estimate of £]( = E jX^l) — . 
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10 ^ 



—©- 

— Euler 

- B- 

— tamed Euler 


— stopped tamed Euler 

— 

- Orderlines 0.2, 0.1, 0.05, 0.01 


Q -2 I..I..I..I..I..I.. 

10° lO'' 10^ 10° lO"^ 10° 10° 

number of time steps 


Figure 1. Error vs. number of time steps 

Figure shows, on a log-log scale, the plots of the error estimates ^ versus the number 

of time-steps n G {2°, 2^, 2°,..., 2^®, 2^°}. Additionally, the powers n“°'°^, n“°'°°, n“°'^, n“°'° 
are plotted versus n G {2°, 2^, 2°,..., 2^®, 2^°}. The results provide some numerical evidence 
for the theoretical hndings in Corollary that is, none of the three schemes converges with a 
positive polynomial strong order of convergence to the solution at the hnal time. 
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